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Abstract 

An effective Chern-Simons theory for the quantum Hall states with edges 

is studied by treating the edge and bulk properties in a unified fashion. An 

exact steady-state solution is obtained for a half-plane geometry using the 

Wiener-Hopf method. For a Hall bar with finite width, it is proved that the 

charge and current distributions do not have a diverging singularity. It is 

shown that there exists only a single mode even for the hierarchical states, 

and the mode is not localized exponentially near the edges. Thus this result 

differs from the edge picture in which electrons are treated as strictly one 

dimensional chiral Luttinger liquids. 
74.10.-d, 73.20.Dx, 73.40.Hm 
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I. INTRODUCTION 



We consider the large scale physics of the fractional quantum Hall (FQH) liquids with 
boundaries. Our analysis is based on the effective Chern-Simons gauge field theory |l -10|. 

In an approach initiated by Wen [I l| (see also Ref. ||), which is different from ours, 
a one-dimensional edge action is added to the original action in order to assure the gauge 
invariance for the Chern-Simons gauge field for restricted geometries. This has been followed 
by a number of authors PJT2|-|H)| who tried to explain the FQH effect based only on the 
Id theory (edge picture) [ P~B| , [TT[] . It has an attractive feature, if correct, that the celebrated 
Tomonaga-Luttinger liquids are realized at the edges of FQH systems ||l8|-p0|. In those 
theories, however, apparently the global properties like distributions of current, charge and 
electromagnetic fields are assumed to be insignificant. Thus one of the main consequences 
is that the charges and currents are localized near the edges. 



Recently Nagaosa and Kohmoto |2T| examined boundary conditions for FQH liquids and 
it was shown that the gauge invariance is not violated if one imposes a physically suitable 
boundary condition. Thus one needs not to modify the original action as was done by Wen. 
In this way, it was possible to study the edge and bulk properties on an equal footing. They 
showed the fractional quantization of Hall conductance by considering both edge and bulk. 

We begin with the composite fermion picture of the FQH effect in Sec. II. Then the 
Chern-Simons effective field theory to describe the composite fermions and the boundary 
condition are introduced in Sec. III. 



Following the line of Ref. [^y] the charge and current distributions of the FQH liquids 
in steady states are discussed in Sec. IV. It is shown that the effect of the interaction can 
not be treated perturbatively from the non-interacting situations. The approach we are 
taking for this correlated electron problem is totally non-pert urvative. The following two 
geometries are considered: (i) half-plane (subsection A) and (ii) Hall bar (subsection B). 
(i) An exact solution is obtained using the Wiener-Hopf method. This method was used, 
in the same geometry, by Thouless J^J to solve the equations of MacDonald, Rice and 
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Brinkman |30| for the non-interacting integer quantum Hall liquid. 

(ii) The asymptotic behavior near the edges are obtained using the theorems of the singular 
integral equations ( Hilbert transformation). It is shown that there exists only a single mode 
even for the hierarchial states. This mode is not localized exponentially along the boundaries 
and our results are not consistent with the edge picture of the FQH effect. 

It is also shown that the charge distribution is finite over the sample including the edges 
This is indispensable to obtain a well-defined theory, since the electron number density must 



be positive. If the charge density has singularities ||30|| , one can not avoid negative divergence 



in the electron number density which leads to an ill-defined theory. 

II. COMPOSITE FERMION PICTURE OF FRACTIONAL HALL LIQUIDS 

Filling factor is defined by 

(2.1) 

where N e is the number of electrons and is the number of flux quanta. For v = 1 the 
first Landau level is totally filled and the higher levels are all empty. Let us see why Hall 
liquids with inverse filling factor equal to an odd integer might be special. Recall that in 
this case the background magnetic field contains an odd number of magnetic flux quanta 
per electron. 



The first argument we cite is due to Jain. |2^ , [23[| Imagine an adiabatic process in which 



we somehow move some of the flux quanta so that p units of flux are attached to each 
electron. For p even, the additional Dirac-Aharonov-Bohm phase associated with moving 
one electron around another is e lwp = 1 and so the statistics of the electron is unchanged. 
The electrons are now moving in a reduced magnetic field B c g = B — 2ixpn where n is the 
number density of electrons. (Note that in our convention, the unit of flux is 2tc.) The filling 
factor has been increased to z/ fj, given by v~ s = [B — 2irpn)/27m = v~ x — p. For z/ fj = m 
an integer, we have v~ x = p + m _1 and 
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v = (2.2) 
mp + 1 

Thus, fractional Hall systems with v — m/ (mp + 1) (p even) may be adiabatically changed 
into an integer Hall system with filling factor m, as was also emphasized by Greiter and 



Wilczek |[24|| . Note that this argument gives us more than we had hoped for. The case we 
wanted to understand, with v~ l = an odd integer, is obtained for m = 1. 

In this way a Hall liquid with v = m/ (mp + 1) with m an integer and p an even integer, 
is related to the integer Hall liquid. We may thus want to argue that since the integer Hall 
liquid is incompressible the fractional Hall liquid is also incompressible. 

Another argument, historically earlier than Jain's argument, is due to Zhang, Hansson, 
and Kivelson [|l| . Imagine attaching all of the flux to the electrons. Thus, each electron gets 
attached to it an odd number of flux quanta. By our preceding argument, the electrons with 
the attached flux quanta become bosons. We now have bosons moving in the absence of a 
background magnetic field. Since bosons can condense according to Bose and Einstein, we 
conclude that a Hall liquid with inverse filling factor equal to an odd integer is energetically 
favored. 

III. EFFECTIVE FIELD THEORY 

The situation described in Sec. II above may be effectively represented by the Chern- 
Simons gauge field theory coupled with external electromagnetic potential. The Chern- 
Simon term plays a role of attaching fluxes to electrons. By solving the equation of motion 
and the Maxwell equation consistency, we can obtain the charge density, current and poten- 
tial profiles. 

The effective Chern-Simons gauge Lagrangian density in the dual representation [25-27 

is 

£ = ^E K u e^a Ifl d u a JX - ]T (^-A,e^ x d u a IX + j/Wf) , (3-1) 
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where ai^ is the Cherm-Simons gauge field, the integer- valued symmetric matrix K is written 
K = I + pC, where / is the m x m identity matrix and C is the m x m matrix in which 
every element is unity, namely, 



C 



1 1 
1 1 



1 



The matrix K specifies the coupling among the Chern-Simons gauge fields and is also related 
to the filling factor by 

" = E(^~% 
u 

of FQH liquid @. The Maxwell term (l/^/j^/f {f Iflv = d^a Iv - d v a^ ) in Q is 
explicitly written as (2/g)[c 2 ff xy — ff Qx — fj 0y ], where g is the coupling constant and c 
is the velocity of the Bogoliubov mode 0. The vector potential (fi = 0,x,y) of the 
electromagnetic field is coupled to the /i-th component of the charge current density which 
is given by 



(3.2) 



i=i 



where the contribution from I-th conserved current densities is 



J? = —e^ x d v a Ix . 

Z7T 



(3.3) 



Note that the vector potential for the constant external magnetic field B has been 
already taken into account in the structure of the K matrix, and is not included in A^. 
Similarly and the density Jj are measured from their average values in the following 
discussion. 

The Lagrangian density is integrated over a sample S. On the boundary dS we impose 



E J i n *\ds = 0, 



(3.4) 



a=x,y 



where n = (n x ,n y ) is the unit vector normal to the boundary. This boundary condition 
simply expresses the physical condition that the current can not flow through the boundary 



dS. Since it is a physical requirement, it is obviously invariant with respect to a gauge 
transformation —>■ + d^(pj. A remarkable fact is that the Chern-Simons term in 
Lagrangian density ( |3.1| ) is also gauge invariant after integration over S with the boundary 
condition 



IV. CHARGE AND CURRENT DISTRIBUTION 



We consider a FQH liquid with filling v = m/(mp + 1) where m is an integer and p is 
an even integer. To study steady-state distributions of charges and currents, we need the 
equation of motion derived from the effective action for the Lagrangian density ( |3.1| ). It 
reads 



£ KjjJ^j = ■——e fluX d u A\ + Adjr/9 
j 2n 



(4.1) 



which is expressed in terms of the current density only and is gauge invariant. It is explicitly 
written 



(4.2) 



where J M =* [Jf , ■■, J(^] and q =* [1, •-, 1] are vectors with m components. 

Let us investigate the stational distributions of potential and charge in a Hall bar which 



K Sng^dt -87rc 2 g- l d y 




3 X 




(d y A - d t Ay)q 


STtg^dt K %-Kc 2 g~ l d x 




3 V 


1 


-(d x A - d t A x )q 


Sng^dy Sng- l d x K 




J° 




Bq 



is uniform in the y-direction, then one may put dtJj = dtA^ = in 



From the 



homogeneity of the system in the ^/-direction, impose d y Jj = dyA^ = 0. Thus we have 

KJ y + ^d x 3° = -^-d x A q, (4.3) 
9 2tt 

—d x 3 y + K3° = 0. (4.4) 
9 



Combining these, we have an equation for J°: 



K 2 3°- 



i>7TC > o2 T _ ^ q2 



d 2 J" = -d 2 A q, 
9 



(4.5) 



Here the potential A Q is connected to the charge current J°(x) = J2i Ji( x ) through 



Aq(x) — — £ / dx' In \x — x'\ J (x') 



(4.6) 



where £ is a constant having the dimension of velocity [&]. Once J° is determined, the 



current 3 y is obtained from ( |4.3|) or (4.4). 

To diagonalize these matrix equations, introduce an orthogonal matrix 



U 



l/y/2 
l/v^ 



-1/V2 
-l/v 7 ^ 



••• l/v^ -1/^2 

The matrix C is diagonalized as 



^00 ••■ ^ 



••• 







Then (|4~5| ) becomes decoupled equations 



;i + mpyr(x) - 



57TC 



4m 



#J u (x) = — diAoix 



(J I Jrn) 



iTCC 



filial J; 
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for / = 1, 



m — 1. 



(4.7) 
(4.8) 



The densities J® — (I = 1, • • • ,m — 1) which are orthogonal to J°(x) and sometimes 
called as "neutral modes" are unphysical degrees of freedom, since they do not have electro- 
magnetic couplings. 

Therefore there is only a single mode even for the hierarchical FQH states and it is not 
exponentially localized near the edge as shown below. This result is in sharp contrast with 
the claims of Wen [|T^] and Macdonald [T7J that there exist a number of edge branches in 
the hierarchical FQH liquids. 



By solving ( |4.7|) , all the long-range behavior of the electronic density and current can be 
obtained. We have two characteristic length scales 

A ^ 7lC (4 9) 

g(l + mpY 

and 

2 g{l + mp) 2 ^ ^ 

These two scales should be much larger than the magnetic length scale since our starting 
point is the long-range effective theory described by the Lagrangian density ( |3 . f | ) . The 
length scale Ai appears as the localization length of an edge mode in plf . , This edge mode, 



however, does not exist when the Hall conductance is quantized, i.e. when the longitudinal 
voltage drop is zero. In what follows, we will denote the ratio of the two scales as 

.-©"-(*)■■ 

and study the effects of the parameter rj. 

In the Lagrangian density (|3.1|) , the particle-particle repulsive interaction is represented 
by the Maxwell term (more specifically, by the spatial part c 2 ff xy ). Thus if r] = A1/A2 = 0, 
( |4.6| ) and ( |4.7p represent a non-interacting case. Remarkably these equations with 77 = are 



essentially identical to those of MacDonald, Rice and Brinkman J30J obtained to study the 
charge and potential profiles in the integer quantum Hall effect. It is quite unexpected since 
our method is based on the effective field theory and is completely different from theirs. 

A. Half-plane 

The Wiener-Hopf method is used to obtain the charge, current and potential profiles of 
the FQH liquids of an infinitely wide Hall bar with a single edge located at x — 0. For the 



integer quantum Hall liquid i.e. rj = 0, Thouless [^9|] analytically solved the equations of 
MacDonald et al. by the Wiener-Hopf method in the same geometry. It is shown, however, 
that the nature of the solutions are completely different for the two cases: 77 = and 77 ^ 0. 
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So the effect of the interaction can not be treated perturbatively from the non-interacting 
situations (rj = 0). The approach we are taking for this correlated electron problem is totally 
non-perturbative. 

In the present geometry, ( |4.6| ) and ( |4.7|) are written 



J°(x) - \ld 2 x J\x) = ^-d 2 x A (x) (0 < x < oo), 

/•OO 

Ao(x) = -£ dx'\n\x - x'\J°(x'). 
Jo 



(4.12) 
(4.13) 



To apply the Wiener- Hopf technique, we separate A (x) as A (x) = Aq(x) + A (x), where 



A+(x) 



A~(x) 



Aq(x) for x > 0, 

for x < 0, 

for x > 0, 

An(x) for x < 0. 



(4.14) 



(4.15) 



Extend (4.12|) to the region x < as 



A 5 



J°(x) - A^J°(x) = ^dt(A+(x) +9(-x)(A+(0) + A+'(0)x)) 
= -MAt{x) - 5'(x)A+(0) - 6(x)A+'(0)x) 

7T4 



where 9{x) is the step function: 9{x) = 1 for x > and #(x) 
denotes differentiation with respect to x. 

Introduce Fourier transforms of these functions as 



(4.16) 

'-oo < x < oo), (4.17) 



for x < 0, and prime 



/(*) 
9±{x) 



/ „\Akx 



dxJ (x)e 



-oo 
oo 



dxA±(x)e lkx . 



(4.18) 
(4.19) 



Note f(k) and g+{k) converge in the upper complex k plane and g~(k) converges in the lower 
complex k plane. Regularize the Fourier transformation of the logarithm by multiplying a 
dumping factor e~ a \ x \ as 

-2 



dxe a|a:| In \x\e 



ilk 



k 2 + a 2 



1 k 
aC H — aln(&; 2 + a 2 ) + A;tan _1 — 

2 a 



(4.20) 



where 7 = lim m ^ OC) (X]™ = i ~ — lnm) is Euler's constant. The r.h.s. becomes — m if 7^ in 
the limit a — >■ and it diverges if k = 0. By discarding this divergence, we have regularized 
equations 

/(*) + A? A; 2 /(fc) = ^(-*V(*) - <(0) + iA+(0)k), (4.21) 
+ <7_(A;)=^/(A:). (4.22) 

Elliminating g+(k), we have 

/(A;) (l + A 2 \k\ + \\ k 2 ) = —Ak 2 g^k) - A+'(0) + iA+{0)k). (4.23) 

In order to analyze ( |4.23| ) we consider factorization of the analytic function 

1 + M+^ 2 = ^|, (4-24) 

where \z\ is defined by \z\ = z if Im(z) > and \z\ = —z if Im(^) < 0; N(z) is analytic on the 
complex plane except for {z = iy\y > 0} and D(z) is analytic except for {z = —iy\y > 0}. 
Solutions for D(z) and N(z) are 

D(z) = rT 1/2 exp - / % dy- , (4.25) 

\ 7r Jo iy + z J 

, n /— 1 r 00 y tan -1 - r JL - 2 - \ 

JV(z) = r/ 1 / 2 exp — / z dy- ^ (4.26) 

\ 7r jo iy — z J 

where a branch of tan -1 y is chosen as < tan -1 y < n. On the real axis (z — k £ R), -D(fc) 
can be written 

D(k) = (1 + + rjk 2 )-V 2 eM k *\ (4.27) 

where 

/ . \ fc Z" 00 , tan" 1 ?/ -tt/2 7T 1 . , 

V(k\r l ) = -J Q dy | - + -(fl(fc) - 1/2) for 77 = 0, 

£• /-oo tan -1 : y 2 — 7r 

= - / d V 2 iz T2 + < e ( k ) " V2) for n > 0, (4.28) 

7T Jo y + K 

Note that the phase factor (p(k;r)) has a singularity at 77 = 0. If 77 = 0, </?(&; 77) varies from 
—71 /A to 7r/4. For 77 > 0, however, it varies from —tt/2 to tt/2. 
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By substituting z = X 2 k to ( |4.24| ) we obtain 

where D(X 2 k) converges in the upper half complex fc-plane and N(X 2 k) converges in the 

lower half complex fc-plane. Rewrite (|4.23|) as 

f(k) ^(k 2 g 2 (k)-A'(0)+tA(0)k) 
D(X 2 k) ~ N(X 2 k) 

then both sides of the equation must be well-defined. Therefore they are entire functions 

which is a constant. Thus we have a solution 

f(k)= const. D(X 2 k). (4.29) 

Now we study the asymptotic behavior of J°(x) for 77 > 0. The differential equation for 
77) is 

d(p(k;rj) _ (3+\oga+k 2 | /3_ logq„/c 2 ^ ^ 



where 



dk 2n(a + k 2 - 1) 2n(a^k 2 - 1) ' 

2 V 2 



OL± 



1 - 2?7 ± VI - 4r7' 

r/(=bl =F 4?7 + VI - 4?]) 



Vi -477(1 - 277 ± vi - 477) ' 

From this, we obtain the asymptotic forms of ip(k; 77) as 

cp(k; 77) ~ — — ((3 + \oga + + /3_ loga_) k - -kilogk — 1) H (A; ~ 0), 

Z7T 7T 
7T 1 , /?_ . \ 1 11., 1 n X 

~ - - — — loga+ + — loga_ - + — -(log - - 1) + •• • (k ~ 00), 
2 Z7T \a + a_ J k r/ir k k 

To study the behavior of J° on the edge, we estimate 

/oo rc 
dkf(k)~ / + k + nk^^idk + ak^ogk- 1)) 

/oo ,111 
dfc(l + fc + t/A; 2 )' 1 / 2 ^- + 4 -(log - - 1)) 

~ [ C dk(l- -k)(c 1 k + c 2 k(logk- 1)) (4.31) 
jo 2 

/oo 111 
^r/- 1 / 2 A;- 1 (c' 1 - + c' 2 -(log--l)) 



< OO. 
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Here we used (|4.27|) , ( |4.29| ) and the constants q, c£ are obtained by combining a± and (3± 
suitably. From this, it is expected that the charge density takes a finite value at the edge 
(x — 0). In other words, J(x) has no divergence if rj > 0. Since we have the Maxwell term in 
the effective Lagrangian which partially comes from Coulomb repulsion between electrons, 
it is a natural consequence. Furthermore, we can obtain the fourier transform of f(k), using 
approximations (1 + k + nk 2 )- 1 / 2 ~ ri~ l l 2 k~ l and (p{k;rf) ~ n/2 — 9(—k)n — c 2 /k, as 

/oo roc ^2 

dkf{k)e~ ikx ~ / dkk- 1 sm(kx + — ) 
-oo JO K 

= const. J (2cx 1/2 ) {x > 0), (4.32) 
= (x<0), 

where Jq(x) = 1 — ^- + |j — • • • is the zeroth Bessel function of the first kind. 



From these and the results of Thouless [ 2~9~ | , we have 



J°(x) ~ X- 1/2 7? = 0, 

~ 1 — const, x rj 7^ 0. (4.33) 



This behavior is shown in Fig. 1. It can be shown from ( f4.4j) that 

— 9 x J y + (l + mp)J° = 0. (4.34) 

Thus the asymptotic behavior of the current distribution near the edge is 

J y (x) ~ —x 1 / 2 + const. rj — 0, 

~ —x + const. n 7^ 0. (4.35) 



B. Hall bar 

For a Hall bar with finite width [-L, L], although time- reversal symmetry (T) is broken 
due to the existence of the magnetic field, the system has TP symmetry where P is parity. 
It leads to J°(x) = —J°(—x) and Aq(x) = — A (— x). Since the integral in (|4.6f) has a finite 
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interval, the Fourier transformation method used in the last subsection for the half plane 
does not work well. Thus we use another powerful method of Hilbert transformation. In 
what follows, we rescale x such that the interval [— L, L] becomes [—1, 1]. Thus we should 
rescale the parameters as Xi/L. We will write these rescaled parameters as Aj for simplicity. 
The derivative of ( }4.6| ) is represented by the Hilbert transformation as 

d x A (x) =^.y.J^dy^- = ^T x [J°(y)], (4.36) 

where T x denotes the Hilbert transformation 

%[f(y)\ - P-v. /' (4.37) 
J-i 7T y — x 

Here p. v. / dy denotes the principal value integral. From (|4.6|) , (|4~7| ) and ( |4.36| ) the equations 
for currents J°{x), and J y (x) are 

J°(x) - X\dlJ\x) = X 2 d x T x [J°(y)}, (4.38) 

^jy(x) = -Xld x J°(x) - X 2 T x [J°(y)}. (4.39) 
c 

i) Ai = (?7 = 0): Suppose that the charge distribution J°(x) has singularities at the edges 
x = ±1, as it has in the half-plane Hall bar. The singularities must be integrable, namely, 
J°(x) ~ (1 ± x)~ a with < a < 1. It can be shown from Theorem II in Appendix that if 
a 1/2, the r.h.s. of ( f4.38| ) has singularities ~ (1 ±x) _a_1 which leads to a contradiction. 



Thus the sigularlity is not allowed except a = 1/2. The factor cot(«7r) in Theorem II 
vanishes if a = 1/2 and the above argument against the existence of singularities does not 
hold. 

Proposition I. // a solution of ( \4.38j ) with X± = (r) — 0) has singularities at the edges 
(x = ±1 ), the power of the singularity must be —1/2. 

Note that this singularity —1/2 coincides with the result for the half-plane (see ( |4.33|) and 

m 



The inverse operation of the Hilbert transformation (Theorem III in Appendix), gives a 
systematic expansion of J°(x) for A2 > 1 as 
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J°(x) = J2K n Jn(x), (4.40) 



ra=0 



where j(x)'s satisfy the relation 

< jo{x) < ji{x) < 

At the edges jo(x) has singularities of power —1/2 and j n (±l) = for n > 1. Theorems I 
and III in Appendix give the explicit solutions for jo and j\ as 

= -7T%' (441) 

jr(x) = --VT^log (l^) . (4.42) 

7T V 1 + X/ 

Using (|4.39| ), we have 



Svr / 1/2 



J y (x) = A 2 1 + — - - v 7 !^ + C?(Aa (4.43) 



p(l + mp) V A 2 \7r 

-1a oW = - ( i±^x - ^-VT^ - + O(A-). (4.44) 

7R; V 7TA 2 2A 2 2A 2 / 



The current distribution ( |4.43j ) is symmetric. It is contrasted with the edge picture in which 



the currents flow in the opposite directions at the two edges. The density, current and 
potential profiles thus obtained is plotted in Fig. 2. 

ii) Ai > (rj > 0): In this case the singularities of the charge distribution at the edges is 
suppressed due to the second derivative term in the r.h.s. of (J4.38[) . 



If J°(x) has power singularities ~ (1 ±x)~ a (0 < a < 1) at the edges, the l.h.s. of ( 4.38[) 
has singularities of power —a — 2. On the other hand the r.h.s. of (|4.38|) has singularities 
of power —a — 1. Thus power singulariteis are forbidden. 

If J°(x) has a logarithmic singularities J°(x) ~ log(l — x) — log(l + x) (which we will 
call as a "simple" logarithmic singularity), one can prove that ( f4.38|) does not hold since 



an-n n \ i n , M 2 log(l - x) - log(l + x) 
d x T x [\og[[l -y)- log(l + y)\ = 

7T 1 — X Z 
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For another logarithmic singularities like higher power of logarithm, for example, to find 
explicit Hilbert transformations becomes more cumbersome. However, it is expected that 

logarithmic singularity 
o x T x I logarithmic singularity ~ = , 

1 — x 1 

from a power counting argument. The second derivative of this logarithmic singularity gives 
a singularity with power —2. Then J°(x) can not have any logarithmic singularities which 
is a contradiction. Thus we assert 

Proposition II. A solution of ( \(.38[) with Ai > (rj > 0) has neither a power singularity 
nor a (simple) logarithmic singularity at the edges (x = ±1^. 

Since the Coulomb interaction repels particles each other, a divergent singularity at an edge 
is physically unacceptable. This observation is consistent with Proposition II. Thus we claim 
that the charge density is finite at the edges. Note that there is a singularity at the edge 
in the noninteracting case r] = ( p9|). This is caused by the absence of particle repulsion 



and disappears once an interaction is taken into account. 
In order to obtain J°(x) we expand it in powers of A2 as 



n=0 



By substituting this into ( |4.38|) , we obtain 



. sinh(x/Ai) . 

smh(l/Ai) 

and 

h(x) - \\d 2 Mx) = d x %\jo(y)}. (4.46) 
The solution of the inhomogenious differential equation (|4.46|) is 

1 



8vrAf sinh(l/Ai 
sinh(x/A 



x 



sinh(l/A 



4 ((-2 - A ie - 2 / Al + 2A 1 )e- 1 / Al Ei f f- 
1) I VAi 
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+(2-A ie 2/Al +2Ai)e 1/Al Ei — 

V Ax 

+2(-Ai cosh(l/Ai) + 2 sinh(l/A 1 ))Chi (J^j 



-2(-Aisinh(l/Ai) +2cosh(l/Ai))Shi ( — 

,Ai 



+ {(-2 - A ie - 2 / Al + 2A X ) (V Al Ei ~ e'^Ei (±± 



+ (2 - A ie 2 / Al + 2Ai) ( e-^Ei ( — ^ ) - e x/Xl Ei 



Ai 



. fl + x 

AT 

1 — X 



(4.47) 



Ai 



+2 <^ (-Ai cosh(a;/Ai) + 2x sinh(>/Ai)) Chi 



1 — x 



Chi 



.fl + x 



Ai 



Ai 

+(-A 1 sinh(x/A 1 ) + 2xcosh(x/A!)) (Shi f -^-^ + Shi / J 



, (4.48) 
Ai / \ Ai / / J J 

where Ei(x), Shi(x) and Chi(x) are the exponential, hyperbolic sine and hyperbolic cosine 

integrals, respectively, which are explicitly written 



Ei(x) 
Shi(x) 



oo e ~t 



"P- v - / -r dt ' 

-X t 



x sinh t 



dt. 



h t 

. / \ i r cosh * - 1 7, 

Chi(xj = 7 + In x + / at, 

jo £ 



(4.49) 
(4.50) 
(4.51) 



and we have used the Hilbert transformation 

sinh(?// Ai) 

1 — x 



sinh(l/Ai 
sinh(a;/ Ai) 



Chi 



7r sinh(l/Ai 
cosh(x/A^ /. 



7rsinh(l/A ly i \ 
Using fl4~6j) and (|4.39|) , we have 



Ai 
1 — x 



•Chi 
Shi 



. /I + x 



Ai 
1 + 2 



(4.52) 



57T 



g(l + mp) 
1 , 

I 



jv(x) = -\\d x (j (x) + X 2 jt{x)) - \2%[Uy)} + 0(X 



2ji 



-p.v. J dx'log \x' - x\ (jo(x') + X 2 ji(x')) + 0{X\). 



(4.53) 
(4.54) 



Note that the current density J y in (|4.53|) has no divergence at the edges and it is again 
symmetric. It is contrasted with the edge picture in which the currents flow in the opposite 
directions at the two edges. The lowest order approximation of Aq(x) is obtained by using 
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-V.v.Jdylog\y-x\j (y) 



Ai 



sinh(l/Ai) 

x j-cosh^/Ax) (Chi (j-^j - Chi (iy^ 
— sinh(a;/Ai) (shi ( — - — \ — Shi ^ 



+ cosh(l/Ai) log ( - — - ) \ . (4.55) 



Ax J V Ax 
- x s 

The electron density is plotted in Fig. 3 (a) in the lowest order approximation (A2 = 0) 
( |4.45| ) and (b) with the correction (|4.48|) . The current distribution is plotted in Fig. 4 (a) 
in the lowest order approximation (A2 = 0) and (b) with the correction (|4.53|) . 

Note that in the lowest order approximation, both the electron density and the current is 
localized near the edges exponentially with the localization length Ai This behavior, however, 
completely disappears again once the corrections are taken into account as seen from Fig. 3 
(b) and Fig. 4 (b) where the bulk currents do not vanish. 

The potential is plotted in Fig. 5. Note that in the lowest order approximation, Aq(x) 
has a maximum and a minimum near the edges in contrast to the non-interacting case 77 = 
( |4.44| ). These maximum and minimum of Aq(x) survive even if the next leading contribution 
is taken into account as seen from Fig. 5 (b). 
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APPENDIX A: THEOREMS OF THE SINGULAR INTEGRAL EQUATIONS 



The following Theorems I, II and III |3T|] are used in the main text. 
Theorem I. 
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i-y 



cot(a7r) 



1 — x 



1 + x) sin(a!7r) 



(Al) 



for — 1 < a < 1. 



Theorem II. Let f(x) be an L p -function (p > 1) which in a small neighborhood (— 1, —1+ 
5) (5 > 0) of the point x = — 1 can be written in the form 

f(x) = A(l + x)- a +g(x) (0<a<l), 

where A is a constant, g(x) vanishes at x = — 1 and satisfies (uniformly) a Lipschitz condi- 
tion of positive order e, i. e. 

\g(x) - g(x )\ < K\x - x \ e . 
Then the Hilbert transform of f(x) has the asymptotic representation 

%[f(y)} = Acot(a7r)(l + x)~ a + 0(1) (x - -1), 
if < a < 1, and the asymptotic representation 

T x [f(y)] = --\og(l + x)+0(l) 

71 

if a = . 

If the point x = — 1 is replaced by the point x = +1, then all remains the same except 
that cot(«7r) is changed into — cot(o;7r) and — log(l + x) is changed into + log(l — x). 



Theorem III. // a given function f(x) belongs to the class L^/s+ e for sufficiently small 
e > 0, the equation 



f(x) = % 



x\y\y )\i 



(A2) 



has the solution 

C 

Vl-x 

where C is an arbitrary constant. 



1 dy l-S- S(y) 



x z y — x 



(A3) 
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FIGURES 



1 The electron density J°(x) for the half-plane, (a) Non- interacting case 77 = in the 
log-log scale, (b) 77 = 0.1, 0.5, 1.0, 1.5 and 2.0. 

2 The density profiles for the Hall bar in the non-interacting case Ai = 0, A 2 = 1.2. (a) 
Electron density J°(x) which has singularities with power —1/2 at the edges x — ±1. 
(b) Current g ^ mp ) J y (x). Note that it is symmetric, (c) Potential ^A (x). Note that 
it is a monotonic function and the electric field is in the same direction in the sample. 

3 The electron density J°(x) for the Hall bar. (a) Ai = 0.1 and A 2 = 0. (b) Ai = 0.1 and 
Ai = 0.15. 

4 The current density g ^ mp ^ J y (x) for the Hall bar. (a) Ai = 0.1 and A 2 = 0. (b) 
Ax = 0.1 and Ax = 0.15. 

5 The potential jA (x) for the Hall bar. (a) Ai = 0.1 and A 2 = 0, (b) Ai = 0.1 and 
Ai = 0.15. Note that they have a maximum and aminimum near the edges. Thus the 
electric field is not in the same direction in the sample and vanishes at these points. 
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